Electric-Field-control of spin rotation in bilayer graphene 
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The manipulation of the electron spin degree of freedom is at the core of the spintronics paradigm, 
which offers the perspective of reduced power consumption, enabled by the decoupling of information 
processing from net charge transfer. Spintronics also offers the possibility of devising hybrid devices 
able to perform logic, communication, and storage operations. Graphene, with its potentially long 
spin-coherence length, is a promising material for spin-encoded information transport. However, 
the small spin-orbit interaction is also a limitation for the design of conventional devices based on 
the canonical Datta-Das spin-FET. An alternative solution can be found in magnetic doping of 
graphene, or, as discussed in the present work, in exploiting the proximity effect between graphene 
and Ferromagnetic Oxides (FOs) . Graphene in proximity to FO experiences an exchange proximity 
interaction (EPI), that acts as an effective Zeeman field for electrons in graphene, inducing a spin 
precession around the magnetization axis of the FO. Here we show that in an appropriately designed 
double-gate field-effect transistor, with a bilayer graphene channel and FO used as a gate dielectric, 
spin-precession of carriers can be turned ON and OFF with the application of a differential voltage 
to the gates. This feature is directly probed in the spin-resolved conductance of the bilayer. 



Graphene has attracted much attention since its first 
experimental fabrication!, due to its exceptional elec- 
tronic properties linked to the Dirac physics of its low- 
energy quasiparticles 2 . Thanks to its extremely high mo- 
bility^, graphene is also a promising material for nano- 
electronics, where however the presence of a semiconduct- 
ing gap is required. One-dimensional graphene-related 
structures like nanoribbons 4 and carbon nanotubes^ have 
been employed succesfully for nanoelectronic devices, but 
their practical applications are limited by the need of 
single-atom precision in the definition of their transver- 
sal width and radius, respectively. Carbon-based materi- 
als like epitaxial graphene on SiC 6 and bilayer graphene 7 
have been shown to be promising for the realization of 
tunneling Field- Effect Transistors (FETs), while the gap 
is not sufficient for conventional FETs 8,9 . 

Spin-orbit coupling plays a crucial role in spintronics, 
providing a way to manipulate electron spin by means of 
an external field. This is at the heart of most proposed 
spintronic devices, such as the Datta-Das spin-FET 10 . 
However, theoretical studies have shown that spin-orbit 
coupling in graphene is extremely smal l 11 ' 12 . Therefore, 
conventional spintronics mechanisms are not applicable 
to graphene. On the other hand, graphene is attractive 
for spintronics because of its long spin-coherence time 1 ^. 
Moreover, simulations show that magnetic doping 1 ^ in 
graphene, or edge functionalization in GNR o 15 ' 16 , lead 
to spin-splitted bands and potentially to a semi-metal 
dispersion relation, particularly attractive for spintronic 
applications^. 

An alternative approach can come from the exploita- 
tion of the interfacial proximity with a ferromagnetic ox- 
ide. Indeeed, it has been proposed that a spin split- 
ting (acting as an effective Zeeman field) can arise in 
graphene due to Exchange Proximity Interaction (EPI) 
between electrons in graphene and localized electrons in 



a FO layer adjacent to the graphene laye r 18 ' 19 . The ef- 
fective Zeeman splitting (which has been estimated to be 
of the order of 5 meV 19 ) acts on the spin of graphene 
carriers inducing a precession around the magnetization 
axis of the FO. The same mechanism was proposed for 
bilayer graphene, and the modification of the electronic 
structure and of the magnetoresistance as a function of 
relative angle between the magnetization axes of the up- 
per and lower FO spacers has been investigated 20 . The 
spin filtering properties of bilayer graphene with multi- 
ple magnetic barriers with EPI were also investigated in 
Ref.l2ll However, the control on the spin-transport prop- 
erties exerted by an electric field, driven by gates, has not 
yet been the subject of investigation. 




FIG. 1: (color online) Picture representing a double-gated 
graphene bilayer structure in which the central part (C) is 
characterized by the use of a ferromagnetic oxide as spacer 
between the upper layer and the top gate. The oxide thickness 
t ox and the interlayer distance do are indicated. The potential 
of the top gate Vgt and the one of the back gate Vqb are 
externally fixed, inducing potential values Vu and Vl on the 
upper and lower graphene layer, respectively, via non-linear 
Poisson equation. 



It is possible to induce an energy gap in bilayer 
graphene by applying an electric field perpendicular to 
the graphene plane^r— . For small kinetic energy, first 
valence and conduction band wavefunctions are driven 
towards different planes of the bilayer by the applied field. 
This means that electrons in the first valence and con- 
duction band have different probabilities to be found on 
the upper and lower plane. Therefore, we can devise a 
way to tune the spin-rotation of carriers in the bilayer 
graphene by reversing the gate voltage. 

We consider a double gate bilayer graphene FET where 
a FO is used as the insulating layer between the bilayer 
graphene channel and the top gate, while an ordinary ox- 
ide (00) is used as insulating layer for the bottom gate, 
as shown in Fig. [TJ EPI interaction will mainly affect 
electrons on the upper graphene layer. By applying a di- 
rect or inverse differential voltage between the gates, we 
determine whether conduction electrons do or do not feel 
the EPI interaction, and whether the associated wave- 
functions are quasi-localized either on the upper or on 
the lower plane. Consequently, we are able to switch on 
or off spin precession. 
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FIG. 2: (color online) (a) Bilayer dispersion curves (solid 
lines) compared with the Dirac dispersion curves for the up- 
per and lower graphene planes without inter-plane mixing in- 
teraction (dashed line), for Vu = 0.05 eV and Vl = —0.05 eV 
(A = 0.1 eV). (b) Projection of the conduction band states 
onto the upper plane as a function of fe, for increasing A. The 
black solid line corresponds to the case in (a), (c) Dispersion 
curve for a bilayer system with EPI interaction on the upper 
plane with Ez = 20 meV. 
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I. MODEL 



We discuss here in very generic terms the difference 
between an insulating material made of an ordinary ox- 
ide and of a FO. In both materials electrons reside in 
similarly localized wavefunctions, as it is proper of an 
insulating material, but in a FO they will be also char- 
acterized by a majority spin component. If we place 
a graphene sheet in proximity to a FO, rather than to 
a OO, in general we can expect a similar contribution 
for the direct Coulomb interaction between electrons in 
graphene and in the oxide, but a completely different 
contribution from the exchange interaction. Indeed, the 
exchange interaction requires "exchanged" electrons to 
have the same spin orientation, and therefore graphene 
electrons will feel a very different effective exchange prox- 
imity interaction (EPI) for majority and minority spin 
components with respect to the FO. Moreover, while the 
direct Coulomb interaction is long-ranged, EPI requires 
an overlap of the wavefunctions of "exchanged" particles. 
For this reason EPI interaction, as pointed out in Ref. |2Q . 
is essentially limited to the graphene layer placed in di- 
rect proximity to the FO, and it is negligible on more 
distant layers. 

We assume here the simplest situation, in which a thin 
FO layer is deposited between the upper graphene plane 
and the top gate, with magnetization M , and an OO layer 
is instead used as insulator between the lower graphene 
layer and the back gate (Fig. [1]). The electronic states of 
bilayer graphene can be described, near the K point, by 



H=H + H m = V I- 
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where Vq = —q{Vjj 



(1) 
V L )/2 and A = -q(Vu - V L ), 

with Vu and Vl the upper and lower layer potential, 
respectively, and q is the value of the absolute elemen- 
tary charge. In the future, let us denote with the in- 
dex U the variables relating to the upper (U) layer, and 
with the index L those relating to the lower (L) layer. 
% = vp(p x + ip y ) is the kinetic energy operator (with 
p y — > —p y for the K' valley), H m is an effective energy 
term due to the EPI with the ferromagnetic insulators. 
We use the parameters t±_ = 0.39 e V 27 ' 28 , i.e. the bi- 
layer interplane coupling, and v/ f=s 10 6 m/s 2 ^. Other 
interlawr coupling terms are neglected, in the spirit of 
Refs. Il9lj20|j29l . as they would not change the qualitative 
features of the phenomenon described in this work. The 
Hamiltonian acts on wavefunctions of the form 



(2) 



where A, B refer to the two inequivalent carbon atoms on 
the upper graphene layer, A' ,B' to that of the lower layer. 
L x and L y are the channel dimensions along X and Y 
directions. Now we distinguish the two spin components 
along the Z axis, perpendicular to the plane, therefore 
Xx, with X — A, B, A', B', has to be regarded as a 




two-component spinor 



Xx 



_ ( (f>xt 



(3) 



The elements Vjj,Vl, it, t± are diagonal in the spinor 
space, while off-diagonal terms can be due to the pres- 
ence of an effective Zeeman field. If we imagine to put 
the upper graphene layer in contact with a ferromagnetic 
insulator having a polarization on the XY plane, the ex- 
change interaction gives rise to an off-diagonal coupling 
in the spinor space of the kind 



h„, — aMS = 



E z (0 



m x + im v 



(4) 



where M 



n x ,m y ) is the versor of the magnetization 
vector M. For simplicity, we assume for the upper layer 
a similar effective Zeeman coupling h m for the A and B 
sites, while EPI vanishes on the lower plane sites A' and 
B'. 



II. BILAYER GRAPHENE 

The idea of controlling spin rotation of bilayer 
graphene is essentially based on the plane-localization 
properties of the bilayer spinors, when a vertical field is 
applied. In particular, let us assume for a moment no 
EPI interaction, i.e. a — (H = Ho). In this case the 
eigenvalues near the K, or K', point, are given by the 
formula 26 



(E ri ) ± = V Q ± 



\ 




h 2 v 2 k 2 , (5) 



with e\ = 4r + h 2 v 2 f k 2 + l -f, and with rj = ±1 for the 
first (-q = —1) and second (77 = +1) conduction (+) 
or valence (— ) band. The bilayer spinors, for a given k 
and energy E, are obtained by solving the linear system 
(Hq — E)i& = 0. In Fig. [2ja) we show the bilayer disper- 
sion curve, compared with the graphene dispersion curves 
E+ = ±Hvfk + — and E ± = ±hvfk ~ — , obtained by 
Eq. [5] by decoupling the two layers (t± = 0). 

In Fig. Hfb) we plot the projection of the first bilayer 
conduction band states (Pu — \xa\ 2 + lxs| 2 ) on the U 
plane. The behavior of the projection on the U plane can 
be easily understood from the bilayer dispersion curve. In 
fact, at k = 0, states of the first bilayer conduction band 
stand on the Dirac point of the U graphene layer. At 
k = 0, the A and B sublattices are not coupled by the 
kinetic term tt in the Hamiltonian in Eq. [TJ therefore B 
and A' sublattices are not mixed by t± and retain their 
original character. Correspondingly, the spinor of the 
first conduction band will have a 100% weight on the B 
sublattice. With increasing k, the bilayer spinors have 
mixed contributions from the two planes and eventually, 



at sufficiently large k, Pjj tends to a constant value. An 
explanation for this comes from the fact that for large k, 
the bilayer conduction bands essentially originate from 
the mixing of only the n-type part of the Dirac cones for 
the U and L graphene sheets (shown in Fig. Uta)), while 
contributions from the p-like part may be neglected. This 
leads to a two-level system with a fixed energy separation 
of A, as plotted in Fig. [U[a), and fixed tunnel-coupling 
%. The solutions of the two level system are E± = 



Hvfk ±yj + fi which correspond to the asymptotic 
behavior of bilayer conduction bands for large k, and 



Pn = 



1 



1+ £■ 



(6) 



for the first conduction band (77 = —1), which explains 
the plateau in Fig. |5fb) at large k. We note that, for 
a given k, P\j of the first conduction band corresponds 
to Pl of the first valence band, and by reversing the 
potential of both layers one would perfectly exchange the 
projection properties of the two bands. 

In Fig. [2][c), we plot the first conduction and va- 
lence band of a bilayer graphene subjected to EPI in- 
teraction as described by the Hamiltonian Eq. [TJ with 
Ez = 20 meV. When the EPI interaction is taken into 
account, electronic wavefuctions traveling on the U plane 
are subject to an effective Zeeman interaction, basically 
proportional to Pu, that results in a spin splitting of 
the bilayer bands by PjjEz- This proportionality clearly 
appears in Fig. [He) , where we plot the conduction and 
valence bands of a bilayer graphene with a FO deposited 
on top of the U layer, considering a potential energy dif- 
ference of A = 0.1 eV between graphene layers. If we 
reverse A, which can be realized by inverting the bias of 
top and back gates, the spin splitting of conduction and 
valence bands is inverted, as well as Pu o Pl, and the 
Zeeman splitting at small k will vanish in the conduction 
band. In a regime in which small k states are responsible 
for transport through a FO-contacted bilayer region, we 
will therefore have a degree of control over the electron 
spin rotation induced by the effective Zeeman field. 



III. TRANSMISSION THROUGH A FO GATED 
REGION 

To calculate the spin rotation properties of our system 
we analyze the transmission through a region of the bi- 
layer in which the EPI coupling is active, as is the case 
in Fig. [TJ For simplicity, we imagine abrupt boundary 
conditions such that the contact with the FO is limited 
to the upper plane of the graphene bilayer from x = to 
x = L c . We consider an incoming conduction band elec- 
tron from the left side (LS) (x < 0) of given k and there- 
fore energy £&, with a chosen f spin polarization. Elastic 
transmission through the active EPI zone conserves k y , 
because of the space homogeneity along the Y axis, but 
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FIG. 3: (color online) (a) Spin differential transmission 
T s = T-f-f — Tfj, through the C region of Lc = 150ac-c with 
EPI and A = 0.1 eV, as a function of the wavevector of the 
incoming particle, (b) T s for a reverse potential energy dif- 
ference of A = —0.1 eV. 



not the spin, and leads to reflected and transmitted com- 
ponents to the left and right side (RS) respectively (of t 
and J, spin character). 

In particular, the LS and RS are described by the 
Hamiltonian Eq. Q] with Ez = 0. We can here find, 
disregarding the spin which is here conserved, four pos- 
sible values of the wavevector k x compatible with k y and 
energy E: k x and — k x , which are propagating modes, k x 
and — k x , which can correspond alternatively to propa- 
gating modes or to evanescent modes, with a finite imag- 
inary part^. The number of propagating modes (with 
real wavevector) corresponds to the number of intersec- 
tion points of the conduction bands with the horizontal 
line E = Ek on Fig. Uta). The remaining modes are 
evanescent. Therefore the total wavefunctions on the LS 
and RS can be written as 



*l(0) = 



l L r + C IN s, 
'fit 



(7) 
(8) 



where the matrices Cl, Cr and C/at are built from the 
spinor set of the bilayer system in Eq. Q] without EPI. 
s is the spin-polarization of the incoming particle and 
is a vector describing the component up and down with 
respect to the Z axis. Our calculation starts from fully 



polarized incoming particles, for which s = (1,0) T . For 
all matrices, rows run over the sublattices {i — A f, A 1 
, B f, B 4-, A' t, A' 4-, B' t, B' I}, while columns run over 
the left region output modes {j — —k x f, — k x i,—k x f 
, — k x 4.} for Cl, right region output modes {j = k x t 
, k x 4-, k x t, k x 4-} for Cfi, and incoming modes {j = fe x t 
, k x 4-} for Cfjv- The output coefficients are collected in 



r = (n,r_i,fi,f_i) T 



(9) 
(10) 



with the tunneling coefficients ti,t-i,ii,i-i for allowed 
modes in up or down spin orientation, and the reflection 
coefficients defined in a similar manner as r\ , r_i , f 1 , r _i . 
In the central part of the system we have a mixing of 
spin components induced by the effective Zeeman split- 
ting. Solving the secular equation for the Hamiltonian in 
Eq. [I] for a given energy E and in-plane momentum k v , 
we obtain 8 solutions for the wavevector k x = a n with 
n = 1, 2 ... 8. The corresponding modes are described by 
the spinor 'F^ where the M index is used to specify that 
these states are for the system with Zeeman interaction. 
The scattering state in the central part of the system can 
be generally expressed as 



*<7 = 






OntfS* - C c a. 



(11) 



The Dirac equation requires the continuity of spinors 
at the boundary x — and x = L, which is now expressed 
by the following linear relations 



l c a = C L f+ C/jvs, 
C R t= C c Pa, 



(12) 
(13) 



with P describing the phase accumulation of the different 
components of the scattering state by traveling through 
the C region. After elimination of a, the problem is 
reduced to the solution of a linear system of the kind 

Mx = S, with x = I _ I , which can be easily solved 

by standard numerical techniques. In practice, a source 
term S, describing the incoming particle, pumps the lin- 
ear system described through the dynamical matrix M, 
which carries all the information about the transmission 
through the central region, and determines the output 
steady state described by x. 

We consider the transmission of our system, which is 
given by the sum of the outgoing propagating compo- 
nents in the RS. We choose FO with magnetization along 
Y, so that in the C region we have spin-splitted bands (Y- 
SSB), eigenstates of S yi while we inject and detect in the 
LS and RS electron spin-polarized along Z. Note that 
as explained before, with direct gate bias, the conduc- 
tion bands will be spin-splitted, while for a reverse bias 
the spin splitting is negligible. In Fig. [3l we show the 
spin differential transmission T s = T^ — T^i through the 
central region with FO deposited on the U layer, with 
Ez = 20 meV and Lc — 150 ac-c- T s is calculated 



as a function of the wavevector of the incoming parti- 
cle in the LS. Fig. (3Ja) is calculated with a direct po- 
tential energy difference between the graphene planes of 
A = 0.1 eV. Indeed a marked resonance is observed with 
negative values of T s . Electrons of this wavevector are 
transmitted through the barrier with a spin rotation of 
7T. In Fig. EJb), where A = —0.1 eV is reversed, such 
a feature is absent and electrons preferentially preserve 
their spin orientation. The spin-transmission properties 
are therefore dramatically affected by changing between 
direct and reverse bias of the T and B gates. In par- 
ticular this resonance falls into the mexican hat region 
of the upper Y-SSB (Fig. [2jc)), where three propagating 
states are active: two from the upper Y-SSB and one 
from the lower one. The resonance condition is given by 
the existence of two propagating modes, one from each 
of the two Y-SSBs, for which Ak x L c — it + 2nir, with 
11 = 0, 1, .. (here n = applies). In fact, an incoming par- 
ticle, spin-polarized along Z, can be transmitted in the 
C region as a linear combination of two states from the 
two bands (eigenstates of S y ). These components, trav- 
eling through the C region, acquire a net phase difference 
of 7r, which corresponds to a net spin-flip process. Note 
that the mexican hat-like dispersion makes it possible to 
have two propagating states with large Ak x , allowing the 
fulfillment of the resonance condition with Lc as small 
as 20 nm. Of course, choosing different values for L c 
leads to different positions of the spin-flip transmission 
resonance. 

For an incoming particle of lower energy, only the lower 
Y-SSB contributes propagating components in the C re- 
gion as shown in Fig. Htc) . The overall transmission prob- 
ability T = T-|-|- + T-fi has an upper limit of 0.5, because 
the propagating component is eigenstate of S y , and can 
be seen as a combination of half and half S z spin compo- 
nents. For the same reason, the spin differential transmis- 
sion is close to zero. For an incoming particle of energy 
above the resonance, instead, there is one propagating 
component for each of the Y-SSBs (see Fig.^c)). How- 
ever Ak x between these components is much smaller with 
respect to the resonance case and their phase difference 
accumulated by traveling through the C region is negli- 
gible. Therefore both the spin differential transmission 
and the overall transmission are close to unity. 



IV. CONDUCTANCE 

A readily measurable property of the system is its con- 
ductance. We have therefore calculated the 2D conduc- 
tance of the bilayer system with FO in the C region. In 
particular we are interested in the spin-flipped relative 
conductance X s = Gf^/G, with G = G T | + Gff. This 
is a measure of the efficiency of spin control in the pro- 
posed device. The 2D two-terminal conductance, which 
is defined by dl x = G x dVos, is expressed as 
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FIG. 4: (color online) (a) Spin-flip relative conductance X s , 
calculated at T = 1.8 K, for a system with direct bias of 
A = 0.1 eV, and inverse bias (b). fi refers to the electrochem- 
ical potential with respect to midgap of the bilayer bands, as 
shown in Fig. [2jc). 



with g = 4 (accounting for the valley and spin degen- 
eracy), v x the group velocity along the transport direc- 
tion, f(E) the Fermi-Dirac distribution function and /i 
the electrochemical potential. Integration is performed 
over the Brillouin zone. 

In Fig. |4] we show the spin-flip relative conductance 
X s as a function of the electrochemical potential for 
Lc from 50ac~c to 200ac-Ci at the temperature of 
1.8 K. Fig. BJa) is computed with a direct gate bias 
of A = 0.1 eV, while (b) shows the case with the re- 
verse bias. The resonance present in (a), corresponding 
to an electrochemical potential for which the spin-flip 
affects more than 80% of transmitted electrons, is com- 
pletely absent in (b), where electrons tend to preserve 
their original spin. This calculation clearly demonstrates 
that we are able to control the spin-flip of carriers trav- 
eling through the system by changing the gate bias. 

The total conductance is thermally activated as /x ap- 
proaches the bottom of the conduction band of the LS re- 
gion. As fj. enters the 3-states spectral region in Fig. [5][c) 
the conductance spin properties are dominated by the be- 
havior of the transmission probability in Fig. GHa) . This 
leads, with a direct bias, to a pronounced resonance of 
X s . A fundamental factor for the appearance of this reso- 
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FIG. 5: (color online) Electrochemical potential as a function 
of Vbrn for a double gate graphene bilayer in quasi-equilibrium 
condition, at a temperature T — 1.8 K, calculated for t ox — 1, 
AV G = 0.023 V and t ox = 2 nm, AV G = 0.034 V. In the inset 
the charge accumulated on the U and L graphene layers in 
the two cases is shown. 



nance is that the spin differential transmission resonance 
in Fig.[3ta) is almost isotropic for small k v , similar to the 
LS dispersion curve, deviating only for large k y . There- 
fore, it is possible, with the appropriate electrochemical 
potential, to adjust the Fermi level to this transmission 
resonance. The states relevant to the conductance will 
therefore be quite well collimated on the spin-flip trans- 
mission resonance. As expected, an increase of the tem- 
perature leads to a broader state population, gradually 
blurring away this feature. 



layer state is described by the Fermi-Dirac distribution, 
and the charge it carries can be distributed on the U and 
L plane according to Pjj and Pl ■ The two descriptions of 
the system, electrostatics and statistics, should be con- 
sistent and their simultaneous solution fixes the U and L 
potentials, and therefore //. 

We focus our calculation on a system with A ss 0.1 eV, 
and analyze the control on the electrochemical potential 
with respect to the midgap of the graphene bilayer. In 
Fig. [5j we show /iasa function of the average gate po- 
tential Vcm, where a potential difference Vgb — Vgt — 
0.023 V and 0.034 V has been applied for t ox — 1 and 
2 nm respectively (values which lead to Vl — Vu w0.1V), 
When the device is empty the electrochemical potential 
linearly increases with Vgw As the electrochemical po- 
tential reaches the bottom of the conduction band, we 
can observe an abrupt change of slope. As the charge ac- 
cumulates in the device, the variation of the electrochem- 
ical potential becomes more difficult, due to the increase 
of the quantum capacitance of the system^!. The spin- 
flip resonance region is easily reached with the tight dou- 
ble gate structure adopted here, which optimizes the elec- 
trostatic control. The considered oxide thicknesses are 
obtained with state-of-the-art semiconductor technology, 
and high-dielectric-constant oxides (the so called high-K 
dielectrics) can allow even better electrostatic control. In 
the inset of Fig. [5] we show the charge accumulated on the 
U and L plane. The charge shows an activation behavior 
in correspondence with the value of Vcm for which the 
electrochemical potential reaches the conduction band. 



SELF-CONSISTENT ANALYSIS 



VI. CONCLUSION 



To provide an indication of the real control that the 
gates exert on the system, and therefore of the observ- 
ability of the phenomenon, we performed a self-consistent 
electrostatic analysis. Indeed, we can fix the absolute 
value of the chemical potential, but we cannot set the 
difference between the electrochemical potential and the 
bilayer graphene midgap. In other words, the value of 
the potential of the U and L layers of the graphene bi- 
layer is the result of the self-consistent calculation, which 
depends on the gate voltages, taking into account the 
capacitive coupling with the gates. We study a dou- 
ble gate FET in which we can independently fix the top 
and back gate voltages {Vgt and Vgb)- Alternatively 
we can give the average gate potential Vcim = Vgt + Vgb , 
which is responsible for rigidly shifting the bands (and 
therefore varying the electrochemical potential // with re- 
spect to midgap) , and the gate voltage difference AVg = 
Vgb — Vgt which opens up the semiconducting gap of 
the graphene bilayer. To describe the electrostatics of 
the system we apply to the graphene bilayer the plane 
capacitor model described in Refs. [aKJCl Another way to 
describe the charge on the U and L plane is band filling. 
In fact, the occuppation of each one of the graphene bi- 



Wc have demonstrated that bilayer graphene FETs, in 
which a ferromagnetic insulator is used as a gate dielec- 
tric, is an interesting system for spin manipulation. In 
particular, we have shown that a good electric control of 
spin rotation can be achieved even in a 2D system with- 
out lateral confinement, at low temperature. We show 
that by switching between a direct and a reverse gate 
polarization, we can modulate the ratio of spin-flipped 
transmitted carriers from more than 80% to less than 
20%. Therefore, the system itself acts as a tunable spin- 
flipping device and offers the possibility to devise spin- 
FETs based on bilayer graphene, exploiting the exchange 
proximity interaction with a ferromagnetic insulator, in- 
stead of the rather weak intrinsic spin-orbit coupling. 
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